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Abstract: 

We investigate in detail a model where the curvaton is coupled to the Standard Model 
higgs. Parametric resonance might be expected to cause a fast decay of the curvaton, so 
that it would not have time to build up the curvature perturbation. However, we show 
that this is not the case, and that the resonant decay of the curvaton may be delayed even 
down to electroweak symmetry breaking. This delay is due to the coupling of the higgs 
to the thermal background, which is formed by the Standard Model degrees of freedom 
created from the inflaton decay. We establish the occurrence of the delay by considering the 
curvaton evolution and the structure of the higgs resonances. We then provide analytical 
expressions for the delay time, and for the subsequent resonant production of the higgs, 
which ultimately leads to the curvaton effective decay width. Contrary to expectations, it 
is possible to obtain the observed curvature perturbation for values of the curvaton-higgs 
coupling as large as 10 -1 . Our calculations also apply in the general case of curvaton decay 
into any non Standard Model species coupled to the thermal background. 
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1 Introduction 

In the curvaton scenario, the primordial density perturbations originate from an additional 
scalar field instead of the inflaton [1]. During inflation the curvaton is just a spectator 
field with no influence on inflation. After inflation the curvaton field starts to oscillate in 
its potential and eventually decays. In the process it imprints its inflationary perturbation 
on the decay products. This perturbation is then converted into an adiabatic curvature 
perturbation by thermalisation. 

The magnitude of the final curvature perturbation depends crucially on the time of the 
curvaton decay. In order to obtain the observed perturbation amplitude of £ ~ 10~ 5 , the 
curvaton condensate (oscillating, homogeneous field) must be relatively long lived. Usually, 
the decay is treated in a simple phenomenological manner and the decay width is simply 
tuned such that the correct perturbation amplitude is obtained. In most cases in the liter- 
ature, the decay of the curvaton is either assumed to be perturbative, or parameterised by 
an effective decay width with no concern for its connection with particle physics (see [2] for 
a recent review). We find neither of these approaches entirely satisfactory, as the curvaton 
can always decay non-perturbatively in a similar manner to (p)reheating after inflation. 
This has been shown in a general case [3], with predictions for the non-Gaussianity pa- 
rameters discussed in [4] . Some aspects of curvaton non-perturbative decay have also been 
discussed in [5], and the case of a self-interacting curvaton has been studied numerically 
on the lattice [6]. In the light of forthcoming results from both Planck and the LHC, it is 
essential to understand the dynamics of models which can produce the observed density 
fluctuations. The predictions for non-Gaussianity in particular may be strongly affected 
by the details of the decay process. Hence, it is of great importance to study the curvaton 
decay in dynamical detail. 

Our aims are twofold: (i) to show that when the curvaton decay products are coupled to 
the thermal background, this completely changes the curvaton decay dynamics as originally 
discussed in [3], and (ii) to make a detailed calculation of the model predictions, taking 
into account all relevant physical phenomena. To achieve these aims, it is essential to 
have a well defined model, where the couplings of the curvaton and its decay products 
are known. For this reason, we choose the curvaton to couple directly, and only, to the 
Standard Model (SM) Higgs boson. Thus, the physics is just the SM plus one additional 
singlet scalar field, the curvaton. In what follows, the curvaton-higgs coupling constant is 
kept as a free parameter. 

The main difference compared to calculations of parametric resonance during inflaton 
(p)reheating [7-11] is that the curvaton condensate evolves in the presence of a thermal 
background that consists of the inflaton decay products. We assume that these include 
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the SM degrees of freedom. As we will see, this has a substantial impact on the resonant 
production of the higgs, effectively blocking the curvaton decay until late times. This late 
decay makes it possible for the curvaton to produce the observed curvature perturbation 
amplitude. 

The paper is organised as follows. In section 2 we present our model and discuss 
qualitatively the time scales of various features relevant for the non-perturbative decay 
of the curvaton. This overview of the dynamics should make it easier for the reader to 
follow the calculations in subsequent sections. In section 3 we solve the curvaton equations 
of motion for both a radiation dominated universe and a curvaton-dominated matter-like 
universe. We find it convenient to parameterise the evolution in terms of a dimensionless 
time variable j, which in a certain regime represents the number of curvaton crossings 
around zero. In section 4 we address the issues relevant for the parametric resonance by 
demonstrating that in the absence of thermal corrections, non-perturbative effects would 
make the curvaton decay very rapidly in a broad resonance. Sections 5, 6 and 7 contain 
the main results of the paper — we turn on thermal corrections and show that these block 
the resonance. Depending on the model parameters, the curvaton condensate can become 
very long lived, surviving until electroweak symmetry breaking (EWSB). We treat narrow 
and broad resonance, both in matter and in radiation domination. From there we calculate 
the time scale of energy transfer from the curvaton to the Higgs field. In section 8 we take 
the results of the previous sections and find out under which circumstances the correct 
amplitude for the curvature perturbation is obtained. Finally, in section 9 we conclude. 

2 Overview of dynamics 

We consider a curvaton a with a quadratic potential interacting with the SM higgs as 



The Higgs field is an SU(2)l doublet and can be written as <3? = 4> T + i(f>ja J , with <j) a four 
real degrees of freedom, X the 2x2 identity matrix, and Gj the Pauli matrices. This gives 
$t$ = \Y Ja <f>%, where a = 0, 1, 2, 3. Where necessary, we take v = 246 GeV and A = 0.13, 



GeV [12, 13]. The free parameters in the potential (2.1) are the curvaton mass m a and 
the curvaton- higgs coupling g. These, together with initial curvaton amplitude cr* and the 
Hubble rate at the end of inflaton reheating, denoted as H*, form the set of independent 
parameters characterising the model. 

Despite the simple form of (2.1), the dynamics of the system after inflation are com- 
plicated, with many different outcomes depending on the parameters. Thus, to aid un- 
derstanding of the paper, this section provides an overview of the various processes that 
can occur as the curvaton field evolves and decays. Detailed calculations and discussion of 
results will be given in sections 3-8. 

Note that in this model the only coupling between the higgs and the curvaton is 
<7 2 <7 2 3>t$, which is the only renormalisable coupling of a singlet scalar to the SM. There- 





as suggested by the LHC discovery of a higgs-like boson with mass 
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fore, perturbative decay of the curvaton does not occur at tree level, at least not before 
electroweak symmetry breaking. Thus, we only need to consider non-perturbative effects. 
These occur because after inflation the curvaton is oscillating in its potential with a fre- 
quency given by its effective mass. As a consequence, there will be a resonant production 
of higgs particles out of the energy stored in the curvaton condensate. 

As in the case of inflaton (p)reheating, the resonances occur for certain momentum 
modes and can be labelled either as broad or narrow resonances. The resonance parameter 
is defined by 

*) = (f^Y, (2.2) 

where £(t) is the amplitude of the curvaton condensate at time t, and m its effective 
mass (which can also depend on time for certain model parameters). A broad resonance 
has g > 1 whereas a narrow resonance has q <C 1. The nature and effectiveness of the 
resonance thus depends on the curvaton initial conditions, its mass, its coupling to the 
higgs and the subsequent evolution. 

In the absence of a thermal background, the curvaton would undergo broad resonance, 
and the timescale of energy transfer from the curvaton to the higgs would be fast (see 
section 4.2). However, this is in general no longer true because of thermal effects, as we 
discuss in detail in sections 5, 6 and 7. 

The basic assumption is that by the time curvaton oscillations start, the inflaton has 
decayed into SM degrees of freedom and that these, including the higgs, have thermalised. 
At this stage, the universe is radiation dominated. Therefore the higgs acquires a large 
effective thermal mass which blocks its resonant production, 

m 2 H (t)^g 2 T T 2 (t) , (2.3) 

where g 2 ~ 0.1 [14] is the effective coupling of the higgs to the SM degrees of freedom 
in the thermal bath of temperature T(t). The curvaton also acquires an effective thermal 
mass given by g\J ($t$) ~ gT(t). 

As the curvaton oscillates, its energy density scales as that of a non-relativistic matter 
fluid, and its relative contribution to the total energy density increases. Before it decays, 
it may or may not begin to dominate over the radiation energy (see section 3.2), but in all 
cases the background radiation remains in thermal equilibrium. Meanwhile, the curvaton 
is oscillating and going through many zero crossings. If the higgs was not coupled to the 
thermal bath, then the decay of the curvaton would be a fast broad resonance, and it would 
occur while the curvaton energy density is subdominant: 

3H 2 (t)M 2 P ~ p md (t) > p a (t), (2.4) 

where p a and p Ta d are the curvaton and radiation background energy densities, H(t) the 
hubble rate and Mp ~ 2.44 • 10 18 GeV the reduced Planck mass. However, the thermal 
blocking means that in many cases the non-perturbative resonant decay is kinematically 
forbidden for a long period of time. The thermal corrections also change the curvaton's 
motion, causing it to oscillate faster. At the time when the blocking is lifted, in most 
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cases the initially broad resonance has become a narrow resonance (q <C 1; see section 6), 
although for large g the resonance could still be broad (see section 5.2). It is also possible 
that before the curvaton decays, it has become the dominant energy density, 



and the universe is effectively matter dominated. As a consequence of the thermal block- 
ing, the nature of the parametric resonance is very different from the naive assumption, 
occurring later and often through narrow resonance. Although decay via narrow resonance 
always begins before EWSB, in many cases (for g < 0(1(T 10 )) the transfer of energy 
from the condensate to the thermal bath is inefficient and not complete before EWSB 
(see Fig. (2)). Resonant production of higgs particles and curvaton decay after EWSB are 
beyond the scope of the present paper but will be addressed in future work. 

The nature of the resonance, the expansion history of the universe, and the actual time 
of the non-perturbative curvaton decay all depend on the parameters m a and g, as well as 
on the initial values of the curvaton field <r* and the hubble rate , at the end of inflaton 
reheating. All these parameters determine the relative magnitude of the curvaton and 
radiation energy densities, and thus the curvature perturbation amplitude (see section 8). 
A schematic timeline of the various events is shown in Fig. 1, but the precise ordering of 
the events depends on the model parameters. The notation used in later sections for the 
dimensionless time variable j is shown in Fig. (1) . 

To chart all the possibilities, we present five separate calculations. First is the case of 
broad resonance in radiation dominated background where the higgs is not in a thermal 
bath (section 4). This is done for demonstrative purposes and to introduce the notation 
needed for the treatment of non-perturbative parametric resonance. We then include the 
effect of the thermal bath and find the Higgs field evolution for both the broad and the 
narrow resonance, in both radiation domination and matter domination. This is done in 
sections 5, 6 and 7. 

Although our calculations are performed within the context of a specific higgs-curvaton 
model, they are relevant for any scenario where a curvaton-like spectator field couples to 
some scalar field that has thermalised in the background. The only considerations which 
relate solely to the higgs are those that concern EWSB. 

3 Curvaton zero crossings: m(T) = m a 

In this section, we determine the zero crossings of the oscillating curvaton, i.e. the times 
when a(t) = 0, in both radiation dominated (RD) backgrounds, and curvaton-dominated 
matter-like backgrounds [which for simplicity we now refer to just as matter-dominated 
(MD)]. This information is all needed for the description of the parametric resonance. 

After inflation, the homogeneous curvaton a will start oscillating around the minimum 
of its potential, with a decaying amplitude due to the expansion of the Universe. The 
homogeneous curvaton obeys the Klein-Gordon equation 



m 2 {t)M 2 p ~ p a (t) > p ra d(t) 



(2.5) 



a + 3H(t)& + m 2 (T)a = 0, 
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Figure 1. Schematic representation of the dynamics after inflation. Note that the ordering of the events 
in the timeline depends on the parameter of the model. This particular example shows a narrow resonance 
(q <C 1) in matter domination. Notation used for the time scales in section 5 and section 6 is also illustrated. 



where we introduced the temperature-dependent curvaton mass 

m 2 (T) = ml+g 2 T 2 , (3.2) 

in anticipation of section 7. For clarity of presentation, in sections 3 - 6 we consider the 
case where the curvaton effective mass m(T) is constant, i.e. m(T) ~ m a . The case when 
the thermal correction gT dominates is treated separately in section 7. 

Note that in the curvaton model the field should be light during inflation, giving 
m a <C -ff*. In this paper we adopt a conservative upper bound m CT /i?* < 0.1. 

3.1 Radiation domination 

In a radiation-dominated era, the scale factor as a function of cosmic time t is given by 
a(t) = a* [1 + 2H*(t — t*)] 1 ^ 2 , with a* and H* the scale factor and Hubble rate at t = 
(we take t* = to be the time of the end of inflaton reheating). The Hubble rate is 

H(t) = — — . (3.3) 

w 1 + 2H*t 2t y ' 

When the curvaton effective mass is m(T) = m a , the dominant solution to (3.1) is 

a (t) = 2 1 /4r(5/4)cr, _ ^ ^ + m j {2H ^ ^ (3 _ 4) 
[m a t + m a /(2H*)]± 4 

where Ji(x) is a Bessel function of order | and T(5/4) ~ 0.9064. The prefactors guarantee 

4 

that o"(0) = cr*. For m a t > 2, the large-argument expansion of the Bessel function gives 

a(t) ~ S(t)sin(m ff t + 7r/8); S(t) = 0.860 - g * .. . (3.5) 

(m cr ty i i' 1 

Thus, the curvaton oscillates with frequency m CT , crossing zero (i.e. a = 0) every time 

ry, 4. 77T 15-7T 23-7T 

" l(jb ~ 8 ' 8 ' 8 ' • 

In what follows, we express the evolution of both the universe and the curvaton in- 
terchangeably as functions of a dimensionless time variable j or as functions of the cosmic 
time t. These are related by 

mat = 7T (j - - J c^l irj. (3.6) 
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Since we are considering the case m(T) = m a , j coincides with the number of zero-crossings 
of the curvaton. 

In terms of this new parameter j, the curvaton amplitude in (3.5) is 

E(j)~ 0.364^. (3.7) 

The scale factor during RD as a function of j is 

a(j) = vf(i-l) 1/2 ^vf J ' 1/2 ' (3 - 8) 

where the scale factor is normalised at the first zero crossing {j = 1) to a(l) = 1. This then 
gives a* « sj 7tt/4^ H*/m a . Because m a /H* <C 1, there is significant expansion before 
the first zero crossing, i.e. a(l)/a* S> 1. 

3.2 Matter domination 

The above discussion assumes that the curvaton is energetically subdominant for the whole 
of its evolution. This is not necessarily true. The curvaton energy density behaves as a 
non-relativistic fluid and may at some point start to dominate the energy of the Universe. 

In a radiation dominated universe (again with m(T) = m a ), the curvaton energy 
density averaged over oscillations is 

(pAj)Uc « \mlY?(j) « 0.0664^ (3.9) 

and the thermal bath energy density is 

P W (J) = ^T\j) = 3M|i7 2 (j) = (3.10) 

where g* = 106.75 ~ 0(1O 2 ). Therefore, equality between the curvaton and background 
energy density p„{j m ) = p Ia d(j m ) occurs at 



J EQ = 1-31 (^) »!, (3-H) 



where "EQ" denotes matter-radiation equality. For j < j EQ , the universe is effectively 
radiation dominated (RD); for j > j EQ , the universe is effectively matter dominated (MD). 
Equivalently, this equality occurs at a temperature 



T m « 0.606 



( 



(7* 

\Mp 



\ gl /2 V 'Mpm a 

2 / \ 1/2 



^ M \w P )^\t) T * (3 - 12) 

where T* is the initial temperature at the end of inflaton reheating, i.e. the "reheating 
temperature". For example, if cr* = 10 16 GeV then the curvaton crosses around zero 
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j BO ~ 10 s times before matter-radiation equality is reached. Taking m a = 0.01//*, we find 



T E0 ~ 10 5 T*, a drop of five orders of magnitude. 



'EQ 

r 

-EQ 

For simplicity, we assume a sharp transition from RD to MD at T = T EQ . In MD, 
i.e. for j > j EQ , the Hubble rate is then 

H(j)\ = = H{Jeq) . (3.13) 

The Klein-Gordon equation (3.1) can again be solved analytically in terms of Bessel func- 
tions, where the coefficients are given by matching the RD and MD solutions at the bound- 
ary. The large-argument limit is valid even at i EQ , and gives for t > t EQ 

a(t) = 0.0817 ? J cos [m a (t + t EQ /3) + S] , (3.14) 

\m a \i -+- i EQ /o)) 

where 8 is chosen to match (3.5) at the top of the oscillation nearest j = j EQ . The curvaton 
amplitude in MD can then be written as 



S(i) | ^0,0195 , {Mp/j ^] M 0,0149 ( £-) . (3.15) 

' MD (jVJeq + 1/3) \Mp) (j/j eq ) 

For future reference, the scale factor during the MD era, i.e. j > j EQ , is 

=(3/4) 2 /3( i/iBQ + i /3 ) 2 / 3 a( iBQ ) 



'BQ w( JEQ 

where 



a 

1 MD 

~i^ 2 (i/i BQ ) 2/3 , (3-16) 



i 



7 J V cr * / 

4 Non-perturbative effects in the absence of a thermal coupling 

In this section we present the calculation for the resonant production of particles that are 
not in thermal contact with the background radiation. In this case, there is no blocking 
of the resonance, which occurs quickly. The decay of the curvaton occurs in the RD era 
within a few oscillations, much before reaching i EQ . Therefore, in this case the curvaton 
never dominates the energy budget. This case is unrealistic if the curvaton couples to the 
SM higgs, but we present it for the purposes of giving definitions and setting the stage for 
the more physically relevant case. We use expressions for the curvaton's evolution from 
sections 3.1 and 3.2 with m(T) = m a . Because of this, the dimensionless time variable 
j counts the number of curvaton zero crossings. Note that comoving wavenumbers are 
denoted by k, whereas physical wavenumbers are denoted by K = k/a. 
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4.1 Higgs equation of motion 

The equation of motion of the higgs field components <f) a (x, t) is 

4> a + 3H<j> a + |a Y, <t>l ~ ^ 2 + 9 2 Z 2 (t) sin 2 (m a t + |) - V 2 j Q = , (4.1) 

From now on we will drop the non-linear term because non-linearities are negligible if 
the background temperature is well above the electroweak scale ~ 0(100 GeV). We also 
eliminate the friction term 2>H<p a by a conformal re-definition of the Higgs field, Xa = 
a 3 / 2 (j) a . This introduces terms proportional to H 2 Xa an d (a/a)Xa- However, these terms 
can be neglected because they are subdominant at subhorizon scales (k 2 /a 2 3> H 2 ,a/a). 
Around each zero crossing, i.e. around each fixed time t(j), we linearise the interaction 
term as sin 2 (m a t + vr/8) ~ m 2 (t — t(j)) 2 <C 1. In fourier space the equation of motion of 
the higgs modes Xa(k,t) then reads 

Xa + Ututim - + i - XV ") = 0, (4.2) 



where 



and 



^(j^^A-i^A, (4.3) 



A = ^gm a (T* . (4.4) 



The energy scale A is the most important scale for resonant higgs production because it 
controls the higgs energy density (see (4.17) and (5.27)). In particular, A characterises 
the scale below which infrared modes with comoving momenta k < k cu t(j) evolve non- 
adiabatically around the jth zero crossing: 



kcut(j) = a{j)K cut {j) « 7877 j 1 / 2 K cut (j) ~ j 1 / 8 A. (4.5) 
Modes with k k cu t(j) behave adiabatically. Using H* ~ T 2 /Mp, (4.4) becomes 

/ \ 1/2 / \ 1/2 

A - 1/2 fe) (£) <»> 

demonstrating that the cut-off energy scale A is always less than T* . 
4.2 Energy transfer 

The physical cut-off scale characterises the time scale Ai(j') ~ 1/K cu t(j), during which the 
IR modes k < k cu t(j) evolve non-adiabatically around t(j). Using a natural time variable 
x = K cut (j)t the higgs mode equation (4.2) can be re-written as 

Xa+(^(j)+X 2 )Xa = 0, (4.7) 
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where ' denotes derivatives with respect to x and n(J) is defined by 

A homogeneous scalar field oscillating around a minimum can create particles of all 
species coupled to it, by parametric resonance. This is a non-perturbative effect which 
has been widely studied in the context of (p)reheating after inflation [7-11]. In our case, 
the curvaton can resonantly produce higgs quanta. This can be investigated using (4.7). 
The higgs occupation number after the curvaton has crossed j times around zero, n k (j), is 
given by [10] 

n k {j) = W k (j) + [1 + 2W k (j)} n k (j - 1) 



- 2VW fc (j)[l + W k (j)Wn k (j - 1) [1 + n k (j - 1)] sin 6(j - 1) , (4.9) 

where 9(j — 1) is a random phase at t = t(j — 1), and W k (j) an infrared (IR) window 
function given by 

W k (j) = e-^\ (4.10) 

The stochastic nature of 8(j) is both because the Universe is expanding and be- 
cause all the higgs modes oscillate many times between two curvaton zero crossings, i.e. 
K cu t{j) j m a ^> 1 (here we are implicitly assuming broad resonance). However, the stochas- 
tic term averages to zero because 6 is uniformly distributed (see [10] for details). Thus, 
the iterative expression for the higgs occupation number n k is 

n k {j) + ij = [1 + 2W K (j)} (n k (j - 1) + . (4.11) 

The first time the curvaton crosses around the minimum of its potential, higgs particles are 
produced only via the spontaneous creation term [first term on the right hand side of (4.9)], 
i.e. n/%(l) = e _7rft2 ( 1 ). This guarantees that only modes satisfying k(j') < 1 are excited. 
Soon after the curvaton begins to oscillate, the occupation number for modes k < k cu t(j) 
becomes large, n k (j) ^> 1. Then we can write 

n k (j) « e lo ^ 1+2W ^n k (j - 1) (4.12) 

and expand each IR window function to approximate the iterative equations as 

n k (j) =[l + 2e~™ 2 ^] n k (j - 1) = [l + 2(1 - vr K 2 (j))] n k (j - 1) 

« 3 [1 - (2tt/3)k 2 (j))] n k (j - 1) « Se'^) n k (j - 1) 

valid only for j > 2 and < 1. Substituting the definition of n(j) from (4.8) we find 

« tf-i e -*mV</V\ (4.14) 
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where 

2 j 

F(j) = i+-J2 rlz - ( 4 - 15 ) 

Note that in the second line of (4.14) we dropped the term containing the higgs mass scale 
Xv 2 , i.e. we approximated e irP ^J)^ a2 ^3)^ v /^) ?a 1. This is a good approximation when »<A, 
which is true unless the curvaton- higgs coupling is extremely small, g < (v/m a ) (v/a*). 
For example, with a* = 10 16 GeV and m CT = 10 6 GeV, then the approximation is valid 
unless g < 10~ 17 . In the case of small g where the approximation is not valid, the first line 
of (4.14) should be used in all equations where ra&(j) appears, such as (4.16) and (4.17). 
The total energy transferred into the higgs (after j > 2 oscillations) is therefore 

Ph ^ = 2ir 2 a 3 (j) J dkk2n kU)\J k2 / a2 ti) + m lU) > 

W 2^(7) ^ j ^ 2 ^ F °' )(fc/A) V(W 2 + y S2 0-), (4-16) 

where the factor 4 accounts for the number of higgs degrees of freedom, and we have 
introduced the effective mass squared of the higgs averaged over the jth semi-oscillation of 
the curvaton as m? h (j) = (g 2 'a 2 (j)) osc = ^g 2 T, 2 (j). In the case of a broad resonance, the 
decay products are non-relativistic 1 , so a/ (k/a) 2 + 1/2 g 2 Yi 2 (j) rj gT,(j)/\/2 for k < k cut (j). 
The total energy transferred into the higgs is found using (3.7) and (3.8), giving 

see Appendix A for details. This formula allows us to estimate the number of curvaton 
semi-oscillations needed for an efficient transfer of the curvaton energy into the higgs. Thus, 
the number of zero crossings Aj needed until the higgs energy density equals that of the 
curvaton condensate is given by the solution to 

p h (Aj) = Pa (Aj) ~ 0.066<T 2 (Ajr 3/2 A 4 , (4.18) 

where p a is given by (3.9). For an order of magnitude estimate, from (4.17) and (4.18) we 
see that Aj ~ O(log 10 (5'~ 5 )) + O(log 10 (m cr /cr*)). Thus, for any sensible parameter values 
of the model, Aj ranges from few to few hundred curvaton zero crossings. The transfer of 
energy by broad parametric resonance is very efficient. 

Thus we may conclude that if there is no coupling of the curvaton decay products to 
the thermal bath, then the curvaton would decay very quickly through a non- perturb at ive 
broad resonance. The energy transfer would be fast, taking no more than O(10) — 0(100) 
curvaton oscillations to complete. There would not be enough time for the ratio p a /prad to 
grow sufficiently to produce the observed density perturbation (see section 8). Therefore, 
for any value of g, this form of non-perturbative decay would rule out the curvaton model. 

x We explicitly demonstrate this in Appendix A. 2 for a higgs coupled to the thermal bath. However, 
this also applies in the absence of any thermal coupling, because the thermal corrections are subdominant 
versus the interaction with the curvaton. 
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However, the thermal corrections to the higgs mass have a drastic effect on these 
results, completely changing the conclusions that would be drawn about the model. We 
investigate this in the following sections. 



5 Broad resonance in the presence of a thermal coupling 

In the curvaton scenario there is a thermal background of inflaton decay products, which 
we assume to consist at least of the SM particles, including the higgs. As a consequence, 
the Higgs field will gain an effective thermal mass, proportional 2 to the temperature T [14], 

m ff( T )~S , T T > 

where the effective thermal higgs coupling is g 2 ~ 0.1, summed over all SM degrees of 
freedom [14]. This effective thermal mass of the higgs will have a large impact on the 
resonance, as we will show below. In effect, the resonance will be blocked for an extended 
period, making it possible for the curvaton to decay late enough to give the observed 
curvature perturbation. 

Likewise, the curvaton also acquires an effective thermal mass, (3.2). We postpone the 
discussion of this to section 7; this section only considers effect of the higgs thermal mass. 

5.1 Higgs mode equation 

Maintaining the same variables used in (4.7), the higgs mode equation now including the 
higgs thermal mass is 

rf2 X« , / 2, , 2 2,^ T 2 {j) 



ix2 + ^« + ^O0 j gJ,+^x« = o. ( 5.2) 

This has again been linearised around each zero crossing, given by integer values of j. Here 
T(j) is the temperature of the bath at time t(j). The effective frequency uj(x,j) around 
each zero-crossing of the curvaton, can be read directly from (5.2) as 

w V^') = « 8 C7)+fl?« 2 0')(^j) 2 + (5-3) 

where x, k, k cut (j) and K cut (j) are given by the relevant expressions in section 4. 

For most of the parameter space, after inflaton reheating the initial 3 resonance pa- 
rameter is = g 2 a^/4m 2 S> 1, corresponding to a broad resonance regime. However, as 
will be shown below, the resonance can be blocked, so that with the evolving resonance 
structure the curvaton decay can also occur in a narrow resonance regime with q <C 1. 

We now estimate the value of j when the resonance parameter q = 1 (for a definition, 
see (2.2)), which marks the transition between the broad and the narrow resonance regimes. 
If q = 1 occurs in RD, then 

^=DL = H^H • (5 - 4) 



2 Additional temperature-dependent corrections to the higgs potential and hence to the effective higgs 
mass are subdominant [14]. 

3 The system begins with > 1 after inflation if g > 2(m a / a,). 
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If q = 1 occurs in MD, then 

to-i)\ m -^&-%~w(%)#>. (5,, 

5.2 Time scale for the onset of broad resonance in radiation domination 

Let us consider the broad resonance regime q > 1, i.e. > m a /g. As long as the 

adiabacity condition duj/dx <C u 2 is satisfied, the higgs modes will evolve adibatically 
as Xa(x) oc exp (J* dx'ui(x'j\ . However, this condition might be violated during a brief 
period Ax = K cut (j)At < 1, precisely when the curvaton crosses around zero (see (5.3)). 
When adiabaticity is violated, there is a non-perturbative production of Higgs bosons. 
Adiabaticity is violated during a brief period Ax only if du/dx > to 2 is satisfied, giving 

(see (4.8) and (5.3)). The second inequality is the condition for non-perturbative pro- 
duction. The time interval is Ax < 1 so (Ax 2/3 - Ax 2 ) < 0(0.1). The second term 
proportional to v 2 increases the non-perturbative production by making the right hand 
side of (5.6) more positive. In contrast, the term proportional to T 2 (j) acts in the op- 
posite direction. When the right hand side of (5.6) is positive, the evolution of the higgs 
occupation numbers is given by (4.11), but with the IR window function 



W K (j) = exp -7T + g y(j)-^- . (5.7) 



T 2 (j) 

Non-perturbative effects take place if the argument of this exponential is small. If W(j) <C 
1 , the higgs occupation numbers are not excited during a curvaton zero-crossing but instead 
evolve adiabatically. Thus we find that the exact condition for non-perturbative higgs 
production is 

< k 2 < -k 2 cut (j) + At; 2 a 2 (j) - g 2 T a 2 (j)T 2 (j) , (5.8) 

7T 

which corresponds to the exponent in (5.7) being < 1. Assuming that the number of 
entropic degrees of freedom of the thermal bath do not change, then a(j)T(j) = a*T*. 
Using (3.8), (3.10), and ignoring factors ~ 0(1), we rewrite (5.8) as 

j « ) (^Y\,( « V=mi>4^v (5 . 9) 



m(T) J \ ra a J \m(T) J m a g^l 2 \m(T) 

This inequality determines the time until the non-perturbative production of higgs quanta 
is unblocked. 

The Xv 2 term in (5.9) is relevant only if 

x 2> 2<%T 2 g\ m a M P 
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or equivalently, if 

T(j) < —V\v ~ O(10 3 ) GeV. (5.11) 

In practice, this is never relevant for the broad resonance in RD, because well above this 
temperature, either the Universe has entered a MD phase, or the system has drifted into the 
narrow resonance regime (q < 1). Therefore, the \v 2 term can be ignored. This shows that 
electroweak symmetry breaking is irrelevant, and that the blocking of particle production 
would be the same if the curvaton were coupled to any scalar field in thermal contact with 
the background. 

Once the particle production is unblocked, the system will start producing higgs quanta 
out of the curvaton condensate energy, by parametric resonance (see e.g. [10]). From (5.9) 
we see that the particle production and thus the broad resonance is unblocked when 



3 ^ Jbr 



Ird g 4 g* V °* J 



This expression is one of the main results of this paper. Equivalently, as a function of 
temperature, it reads 

T<T BR \ m ^^(^V^~P. (5-13) 

The analysis of this section is only valid under a number of conditions: (i) the universe is 
in the RD regime; (ii) q J> 1, i.e. j BR | ^ Jeq ^^d j BR | < j q =i\ ; (iii) the curvaton 
thermal corrections are negligible, i.e. gT <C m a and thus m(T) ~ m a . In section 7 we 
discuss the effect of the curvaton's thermal mass. 

Once the resonance has begun, it takes a certain number of curvaton zero crossings, 
Aj, to efficiently transfer the energy from the curvaton condensate to the higgs particles. 
The total time for an efficient decay of the curvaton is then j BR | + Aj . In general, 
Aj ~ O(10) — O(10 2 ), so Aj <C j BR | We discuss the calculation of Aj in section 5.4, 
after calculating the time scale of unblocking in a MD regime. We demonstrate the effect 
of the curvaton's thermal mass on the resonance in 7.2. 

5.3 Time scale for the onset of broad resonance in matter domination 

For large g, broad resonance occurs quickly in RD, before j = j EQ . However for smaller g, 
the resonance is blocked until the universe becomes effectively matter-dominated because 
of the oscillating curvaton. In many cases, if the resonance is blocked until the MD era, 
then it has become a narrow resonance. However, broad resonance in MD could still occur 
for certain large cr* and small m a combinations. We therefore cover this possibility for 
completeness. Let us recall again that we are considering m{T) = m a in this section; the 
resulting unblocking temperature changes when the curvaton's thermal mass is considered, 
see section 7.2. 

The temperature T EQ that defines the onset of matter domination, (3.12), has the same 
functional dependence as the threshold temperature for non-perturbative higgs production 
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in RD, T BR (5.13). Their ratio is given by 



1/2 4 



~ ^fc. (5.14) 

r„J a 2 



Thus, if the curvaton-higgs coupling is g -C gl g 2 ~ 0(1CT 2 ), then matter domination 
begins much earlier than non-perturbative Higgs particle production. To estimate the 
evolution of the higgs occupation numbers in the presence of an expanding MD background, 
we must first study the dynamics in this MD background. 

Analysing the higgs equation of motion as in the previous subsection, we arrive at the 
cut-off scale characteristic of the MD era. We find 

~ (ff)' r ' /2A - (5J5) 

The comoving cut-off scale k cut (j) I = a(j)K cut (j)\ is then given by 
WJ)I MD = 0-825 (^) ' [j/j eq + 1/3] 1/6 A 

~(f)V /6 A (5.16) 

where we used a (j)| MD from (3.17). The scaling with the number of oscillations is stronger 
than in the RD era, but still very mild (oc J 1 / 6 ). Note the extra factors (Mp/er*) 1 / 2 3> 1 
and (Mp/cr*) 7 / 6 3> 1, as compared to the analogous scales in RD, (4.3) and (4.5). These 
reflect the growth of the cut-off scale during the RD period. 

From here we can again calculate the time until non-perturbative effects are unblocked. 
The thermal blocking becomes inefficient when 

k 2 - Xv 2 a 2 (j)+gla 2 (j)T 2 (j) > ^L(i)| MD , (5-17) 

all evaluated in MD when j > j EQ . Using a(j)T(j) = a EQ T EQ and neglecting again the Xv 2 
term, the timescale for lifting the thermal blocking is given by 



(Mp\ 



= ^ 4 L+w (5-18) 

This is another important result of our paper. Compared to the RD stage, the extra factor 
(M P /a*) 1/2 » 1 delays the non-perturbative effects, whereas the power ' speeds them 
up. 

Thus, if J B r| rd < Jeq' then the non-perturbative effects take place during the RD era 
and the time that we have to wait until thermal blocking; becomes inefficient is 7' I . The 

I RD 
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total time required for the curvaton to decay efficiently into the higgs is then j BR | + Aj. 
^ Jbr| rd > Jeq then the non-perturbative effects instead take place during the MD era, 
and the total time is j BR + Aj. Note that in some cases, non-perturbative effects can 



occur very soon after the RD/MD transition, i.e. j BR | ~ j EQ . This is the case when 
(from (5.18)) 




In the next section 5.4, we calculate Aj in the different cases. 
5.4 Particle production in broad resonance 

The broad resonance becomes unblocked when j = j BR . At this point, the higgs occupation 
numbers corresponding to long wavelength modes (k < k cu t(j)) will start growing each 
time the curvaton crosses around zero (integer j). Note that Aj counts curvaton zero 
crossings regardless of whether the curvaton's thermal corrections are initially important 
or not. This is because we always have gT BR <C m a when the resonance is unblocked. 
The subsequent behaviour will depend on whether the created higgs particles thermalise 
between consecutive zero crossings. This is a consideration unique to curvaton decay, 
because in inflaton (p)reheating models there is no thermal background. It is an important 
consideration because non-perturbative effects build up new occupation numbers over the 
previous ones. To quantify this possibility, we need to compare the equilibration time 
required by the higgs to thermalise to the oscillation semiperiod of the curvaton. 

The higgs thermalisation rate at time t(j) is r t h(j) ~ 9 T T(j) [15], and thus the ther- 
malisation time is Ai t h(j) ~ g^ T (j). The curvaton semiperiod is nm^ 1 . Thus, we 
conclude that thermalisation of the decay products occurs if m a -C T(j BR ), but does not 
occur if m a > T(j BR ). 

5.4.1 No thermalisation of decay products: m a S> T(j BR ) 

When j > j BR (either in RD or MD), higgs particles are produced at each curvaton zero- 
crossing. The first burst of higgs particles will be created on top of the initial background 
thermal distribution, so n k (j BR ) = ^fc h (r(j BR )), giving the recursive relation 



where Wk(j BR + 1) describes the production of higgs particles out of the vacuum (like in the 
usual inflaton (p)reheating case) and [1 + 2W k (j BR + l)]n k (j BR ) represents the stimulated 
creation of particles over the thermal ensemble (unique to the higgs production by the 
curvaton). The second term in (5.20) describes both a growth of the occupation numbers 
but also a spectral distortion of the original thermal distribution. The following burst of 
particles gives 



Mj BR + 1) = W k(j B n + !) + [! + 2W fc (j BR + 1)] n k (j BR ), 



(5.20) 



n k (j BR + 2) = W k (j BR + 2) + [1 + 2W k (j BR + 2)] n k (j BR + 1) 

= W k (j BR + 2) + [1 + 2T^(j BR + 2)] [1 + 2I^(j BR + 1)] n k (j BR ) + 
[l + 2W fc 0' B R + 2)] W{j BR +l), 



(5.21) 
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again with terms describing the production of higgs particles out of the vacuum and the 
stimulated growth of the existing higgs distribution. In the same manner, a generalised 
expression for the higgs occupation number i zero-crossings after j = j BH is given by 

MJnp +i) = W k (j NP + i) + [1 + 2W k {j NP + i)} n k (j NP + i - 1) 



H(l + 2W k (j NP +a)) 



.a=2 



W k (j NP +!) + (! + 2W k (j NP + 1)) n k (j NP ) 



+ E II ( X + 2 ^0np + c)) W k (j NP + b) . (5.22) 

6=2 \c=6+l / 

In this generalised expression, the second term (containing indexes b and c) describes 
particles created directly or indirectly out of the vacuum. After a few oscillations, this 
term becomes subdominant to the term containing the index a. The term W k (j NP + 1) 
describes the first creation of particles from the vacuum, whereas (1 + 2W(j NP + 1)) n k (j NP ) 
describes the deformation of the original thermal higgs ensemble. 
The higgs energy density is given by 

Ph{j BR + *) = 2ir 2 a 3 (i) / dk fc2 ™*^Ba + *) V^AK^ + m VS) > ( 5 - 23 ) 

where the scale factor a is evaluated at j = j BP + i, the factor 4 accounts for all the 
higgs degrees of freedom, and represents the effective higgs mass at the moment 

j = j BR + i, including both the curvaton-higgs interaction and thermal contributions. 
When Ph(j BB + Aj) = PctUbr + we wm consider that the energy of the curvaton 

condensate has been efficiently transferred into the higgs. 
The curvaton energy density is given by 

Pa (j m + Aj) = ^£ 2 (j NP + Aj) (5.24) 

In principle, Aj could be computed by solving (5.22) recursively, obtaining p\ t (5.23), and 
comparing this to p a (5.24). However, as we show in Appendix A, approximations for n k (j) 
can be used to give the closed expression 

n{l ™ +Al)a e7— A Wf^TAT) (5 ' 25) 



where the resonant parameter q is also evaluated at j = j BR + Aj, and f(q) is given by 

/(?) = 1 + / 2+ J 4 fi ( 5 - 26 ) 

exp ( 5t g 1 / 4 — 1) 

For a detailed derivation of (5.25), see Appendix A. Here we will list some of the key points 
required for the derivation: 

1. The factor 4 is because the higgs has 4 components. 
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The factor /(g) accounts for the fact that the transfer of energy takes place via the 
building up of the occupation numbers both from the vacuum, and from the already 
existing thermal ensemble. 

2\ Aj-1 



The factor (l 
nential in Aj. 



encodes the growth due to the parametric resonance, expo- 



4. The factor (l + ~) ^ 1 differs from the case of no thermal bath (section 4), which has 
the factor 3 Ajl_1 . Thus the resonance in a thermal bath is less efficient ((1 + 2/e) < 
3). The extra factor 1/e is due to the thermal mass term in the window functions 
(see (A.ll)). 

5. As usual, the energy is diluted due to the expansion of the Universe, as described 
by the factor 1/a (j BR + Aj). This factor is a combination of 1/a 3 corresponding to 
the dilution of non-relativistic particles, and 1/a corresponding to the decrease of the 
effective mass mu{Aj). 

6. Compared to the case with no thermal bath, there is an extra suppression from the 



-3/2 

(1+27^)/ 

for a substantial period. 



factor ( 1 + /^[r, }• 



This is because the parametric resonance has been blocked 



7. The broadness and thus the efficiency of the resonance is controlled by the prefactor 
g 1 / 4 , so larger q gives a more efficient resonance. 

8. The expression is valid for both RD and MD backgrounds. Explicit expressions for 
a(j) and k cut (j Bn + 1) can be found in sections 3, 4 and 5, see (3.8), (3.16), (4.5) and 
(5.16). 

See Appendix A for further clarifications. As an example, for the case when broad resonance 
is allowed to occur during RD, we have 



Ph [J 



+ Aj) « 0.004 /(g) q 



1/4 



(l + |)Ai- 



1 + 



Aj-1 
(e/2+1) 



3/2 



+ 1 



1/2 



+ Aj 



A 



(5.27) 



The resonance can be considered complete when Ph(j BR + Aj) ~ p a (j BR + Aj). This 
gives Aj. In general, the time scale for energy transfer in a broad resonance is Aj <C i BR 
in both RD and MD. Once we have calculated Aj, the time scale of the curvaton decay 
can be estimated as td ec = ^Ubb + ^j) m a 1 - From here the effective decay width is 



e// 



1 

^dec 



(5.28) 



5.4.2 Decay products thermalise: m a <C T(j NP ) 

In the case that the higgs particles thermalise between each curvaton zero crossing, the 
process of energy transfer is altered. Between zero crossings, the additional energy in the 
newly-created higgs particles is distributed between all species in the thermal bath. Thus, 
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each time the curvaton crosses zero, there is a universal injection of energy into the thermal 
bath. To see this, note that the higgs occupation number after the first zero crossing (when 
non-perturbative effects are not blocked anymore) is 

"fc(i B R + 1) = Minn) + [1 + 2n fc (j BR )]I^(j BR + 1), (5.29) 

where n k {j BR ) = n| (T(j BK )) is the initial thermal background distribution. Ignoring the 
thermalisation, after the next curvaton zero crossing we find 

n k (j BR + 2) = n k (j BR + 1) + [1 + 2n k (j BR + l)]W k (j BR + 2). (5.30) 

However, because the higgs particles thermalise between the first and second crossings, 
then n k (j BR + 1) can be replaced by the thermal distribution n k h (j BR + 1). Following this 
logic we obtain the general expression 

MJbr + = 4 h B R + i ~ 1) + [1 + 24 h B R + i - l)]W k (j BR + i). (5.31) 

Every time the curvaton crosses zero, a spectral distortion of the type ~ (1 + 2n t k )W k is 
added to the existing thermal distribution n|. h . After every curvaton zero-crossing, there 
is consequently an injection of energy into the thermal background, given by 

Ap(j BR ^) = 27r 2 a 3(j- +i) / dkk 2 [l+2nt{j BR M-l)]W k {j BR +i)^{k/ay + m 2 h (i) (5.32) 

with mh{i) the higgs effective mass at j = j BR + i defined as in section 5.4.1 (also see 
Appendix A). 

The injection of energy corresponds to an incremental increase in the background 
temperature, T 4 (j BR + i) ->■ T A {j BR + i) + (ATj) 4 , so 

Ap(j BR +z)^^x(AT, l )\ (5.33) 

Assuming that the number of degrees of freedom in the thermal bath do not change, we 
use T oc 1/a to give the energy density of the thermal background at j = j BR + i, 

A*„ + o - + o + W (f^) 4 



1=1 

*f T A ( q 0br) X 1 

30^* ™\ a {j BR +i) 



* ( m V a\j BR +l) 

^vw a4 ^ R ) 



(5.34) 



Using j BR S> 1, and provided that i <C j BR , we approximate a(j BH -+- 1) ~ o(j Ba ) in (5.34), 
giving 



7T 2 



30 

with 



^br+^^z)-^, (5.35) 



G(i) = l + ]T(ATi/T BR ) 4 . (5.36) 



i=i 
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Note that AT/ can be obtained from (5.32) and (5.33). In Appendix A we discuss 
approximations for G(i), by considering the small and large momentum behaviour of both 
the window function W&(j) and the thermal ensemble distribution nj, h (j). We find that the 
integrand [l+2n]j. h ] in (5.32) is constant for k < k cut (j BR ) and exponentially decreases for 
k > k cu t(j BR ). This implies that AT/ depends only weakly on the number I of zero-crossings 
after j . Therefore we can approximate G(i) as 

G(i) ~ 1 + i x (AT BR /T BR ) 4 , (5.37) 

where AT BR is a universal function that depends on the resonance parameter q (evaluated 

at Jbr)> s iven b y 

(AT BR ) 4 = ^ [ dkk 2 [l + 2nt(j Bn )}W k (j BK )J(k/a)* + m^W) 
g*ir a \3 BR ) J 

15 ( 2 \ q^KUj^) 



g^ V e 9 T 7 r 1 i 1 -1/ 4vre 
= 0.044 U + ^ ) ^T B 4 R , (5.38) 

where for the last line we use the fact that at the onset of non-perturbative particle pro- 
duction, <7^,Tg = K^ ut (j BR )/TT. Thus, the dependence on T BR in G(i) drops out, and G(i) 
only depends on the resonance parameter q evaluated at j BR . 

Putting everything together, the thermal bath energy density grows as 



p th (j BR + Aj) « p th (j BR ) 



4 1 



(5.39) 



We obtain the number of curvaton semi-oscillations Aj required for efficient transfer of 
curvaton energy into the thermal bath from Pcr(j BR + Aj) = p th (j BR + Aj). We typically 
find that Aj S> 1, and that Aj decreases for larger q. 



6 Narrow resonance in the presence of a thermal coupling 

In the previous sections, we implicitly assumed that the non-perturbative production of 
higgs particles occurs by broad resonance with q > 1. However because q cx S 2 (j) (2.2), the 
parameter q decreases with time. If the resonance is blocked for a sufficiently long time, 
then when the resonance is unblocked, the curvaton has entered the narrow resonance 
regime. For small values of the coupling g, the resonance is only unblocked in a narrow 
resonance regime with g< 1. Therefore, the calculations in section 5 should be repeated 
for the narrow resonance. 

We remind the reader that both the higgs and the curvaton can gain a large effective 
mass in the thermal background. In this section, we still focus on the effect of the higgs 
thermal mass whilst m{T) = m a for the curvaton. In the case where the curvaton's thermal 
corrections are dominant, m(T) ~ gT, the narrow resonance is blocked in most cases, see 
section 7. 
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6.1 Thermal blocking in narrow resonance 

The two resonance regimes have qualitatively different behaviour. In a broad resonance 
the particle production takes place only during an interval At that is much smaller than 
the curvaton semiperiod Tr/m a . The higgs quanta are generated when the higgs mode 
functions evolve non-adiabatically when the curvaton field is close to zero. The ratio 
At/m" 1 is controlled by the resonance parameter: 

At 



T ~<r 1/4 - (6.1) 



lllr 



For a very broad resonance, m^At <C 1. The particles are also created within a broad 
infrared band of momenta k < k cu t(j), with k C ut(j) some model-dependent cut-off scale 
(see (4.5) for a RD background and (5.16) for a MD background). 

In contrast, a narrow resonance is a continuous process, and excites only modes within 
a very thin shell in momentum space. Like the broad resonance, the narrow resonance 
is also blocked by thermal effects until the temperature drops below a certain threshold; 
however, this temperature threshold is different than for the broad resonance. 

Let us now calculate the temperature threshold for the narrow resonance, again ignor- 
ing the higgs non-linearities and the negative mass term — Xv 2 . The higgs mode equation 
including thermal corrections is 

^2 +5 2 X 2 (i) sin 2 (m a t + | 



+ 3H<p a + [ ^ + 5 2 S 2 (t)sin 2 (m a t+ -) + 5 2 T 2 ) </> a = , (6.2) 



where ' represents derivatives with respect cosmic time t. A narrow resonance is also a 
non-perturbative effect [10] but can be understood as a two-to-two body process where two 
curvaton quanta produce two higgs particles. This is in contrast to the broad resonance, 
which is a collective field theory effect. The production rate in narrow resonance must be 
obtained from a field theoretical analysis of (6.2). Because the curvaton is homogeneous, 
we interpret it as a collection of curvaton quanta with zero momentum k = 0. Conservation 
of energy then implies 

2m a = 2E(k) , (6.3) 
where E(k) is the energy of each higgs particle produced in the process, 

E{k) = ^ + Mt)™ 2 a sin 2 (m a t +^)+ g\T 2 . (6.4) 

The resonant parameter, defined as before by q(t) = g 2 T, 2 (t)/4:m 2 , has been explicitly 
introduced. Particle production requires positive energy solutions for fixed k 2 in (6.3). 
Therefore the following inequality must hold: 

vr\ . 2 



g 2 T T 2 + Aq(t)mi sin 2 [m a t + - J < m 2 . (6.5) 

For the narrow resonance, g < 1 and thus the threshold temperature for the onset of 
non-perturbative effects is 

T NR = ^(l + 0( 9 )) . (6.6) 
9t 
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In contrast to the broad resonance, the narrow resonance threshold temperature is almost 
independent of the model parameters. The narrow resonance will always become unblocked 
before the EWSB, provided that m a 3> 100 GeV. However, in some cases, the transfer of 
energy is so slow that no appreciable energy transfer takes place before EWSB. 

Note that we have presented the narrow resonance results only in terms of temperature. 
Expressions in terms of dimensionless time j and other quantities can easily be derived, 



similarly to section 5. For example, if T EQ 
RD and we find 



< T NR then the narrow resonance takes place in 



0.48^ M P 



1/2 

9* 



(6.7) 



In the opposite case of T EQ > T NR , the narrow resonance begins in MD and we find 



0.5 



4^ 

3/8 
J* , 



(Mp\ 
\m a J 



6.2 Particle production in a narrow resonance 

The study of the particle production in a narrow resonance is usually done by mapping 
the field fluctuations (6.2) onto the Mathieu equation. We show this in Appendix B. The 
outcome is that only modes with momenta k within a very thin shell will be excited. In 
particular, the growing higgs modes (conformally transformed as in section 3.1) can be 
described by [10] 

Xk oc e ^ k ' t)m ^, (6.9) 

with fi(k) the Floquet index given by 



\ 



T 



(6.10) 



(see Appendix B for details). The occupation numbers will grow as re/%(i) 



\Xk\ 



x 



e /j,(k,t)m a t_ As can be seen from the structure of fi(k, t), only modes where y(k/a) 2 + g 2 T 2 
is within the band [1 — 3q/2, 1 — q/2]m a will be excited. 

We show in Appendix B that once the narrow resonance has begun, the growth of the 
higgs energy density can be approximated by 



PhUtm + A i) 



f/2 e n g A jm 4 



5.11 



12^ 2 [e T - 1] 

where q is evaluated at j = j NR . This will give the minimum number of zero crossings 
Aj needed to transfer the energy. Comparing the higgs energy density with the curvaton 
energy density (5.24), we find that Aj is of the order 



log (sV /2 (i NR )) 



(6.12) 



Therefore, the narrow resonance is completed at j 
in many cases A j > j NR . 



j NR + Aj. In general, Aj S> 1, and 
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7 Effect of the curvaton's thermal mass m{T) 



In this section, we discuss the effect of the curvaton's thermal mass on the results presented 
in sections 5 and 6. The curvaton's effective mass is (3.2) 

m 2 (T) = ml + g 2 T 2 . (7.1) 

For simplicity we ignore numerical factors of 0(1) in the thermal correction g 2 T 2 . Note 
that the term g 2 T 2 is only valid when the effective higgs mass is small, i.e. ga <T, In a 
non-negligible fraction of parameter space, the temperature correction g 2 T 2 dominates the 
effective mass for some temperatures, so m(T) ~ gT. This occurs when the curvaton-higgs 
coupling is sufficiently large, g > m a /T*. We first show how the dynamics of the curvaton 
and thus the crossings around zero are affected. We then demonstrate the effect on the 
resonance when the curvaton's effective mass is dominated by g 2 T 2 . We find that the 
narrow resonance is always blocked, and that broad resonance can only occur for a small 
fraction of parameter space, always in RD. 

7.1 Curvaton dynamics when m{T) ~ gT 

In the limit g S> m a /T, the effective mass in the curvaton's equation of motion (3.1) is 
m{T) ~ gT, and the solution in RD becomes 

a{t) = E(t) cos |V 1 + 2H *t ~ l] ) . ( 7 - 2 ) 

where 

nt) - Trrm (7 - 3) 

and where cr(0) = a* has been imposed. Note that this solution is valid independently of 
the size of the ratio gT*/H*. As in the case m{T) ~ m a , it describes a damped oscillator. 
However, now the effective frequency is time- varying. Consequently, the curvaton does not 
cross around zero at regular intervals. Using n to denote the number of zero crossings in 
this regime, we find that these are characterised by the condition 



ir (n-|) | 1 
2 gT* H* 



n = 1,2,3,... (7.4) 



In the limit gT* <C H*, the universe expands significantly during one oscillation of the 
curvaton. Thus, the period of oscillation is not constant, and t n is quadratic in n as 
gT*t n w (ir/2)(n — 1/2) 2 (H*/ gT*). However, in the limit gT* 3> H*, the curvaton oscillates 
many times in one Hubble time. Thus, the period of oscillation is almost constant, and t n 
is approximately linear in n as gT*t n ~ 7r(n — 1/2). 

The solution (7.2) is only valid when gT > m a , which translates in cosmic time to 

t < ttran = \ [{gT* / m a ) 2 - l] H^ 1 , (7.5) 



4 Not to be confused with n used in other sections of the paper. 
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where ttran gives the moment when the effective curvaton mass transitions from domination 
by gT to domination by m a . Equivalently, the transition temperature is 

Ttran = — • (7.6) 

9 

As one can see from (7.5), the initial ratio gT*jm a controls how long the curvaton remains 
in the regime with m(T) ~ gT. If gT* 3> m a , the dynamics of the curvaton are described 
by (7.2) for a long period (measured in terms of the initial Hubble time H* 1 ). However, 
the temperature of the radiation background is decreasing as the universe expands. Thus, 
when the temperature drops to T trans = m a /g at t = ttrans, we must match the solutions 
in the two regimes: (7.2) corresponding to the curvaton dynamics for m(T) ~ gT, and 
(3.4) corresponding to the curvaton dynamics for m(T) ~ m a . For simplicity, we match 
the amplitudes exactly at t = ttrans- This matching always occurs in the RD era. 

Note that although we still define m a t = irj, j no longer counts the curvaton zero- 
crossings, which are no longer regular (see (7.4)). 

7.2 Resonance in the presence of a curvaton thermal mass m(T) gT 

In sections 5 and 6 we derived the conditions necessary to unblock the resonance when 
the curvaton's effective mass is m{T) ~ m a . In this section, we consider the regime where 
the thermal mass dominates, i.e. m(T) ~ gT, and thus where the curvaton dynamics are 
described by (7.2). There are two important considerations: the expression for the curvaton 
amplitude £(t) is altered according to (7.3), and the zero-crossings are described by (7.4) 
instead of by j. If g < m a /T*, then the effective mass is dominated by m a , and the results 
from sections 5 and 6 hold. However if g > m a /T*, then the considerations in this section 
are important. We study the RD regime, because the curvaton cannot become dominant 
if the thermal corrections dominate its effective mass. 

Thus, to determine the cut-off scale for broad resonance in RD, we first linearise the 
interaction term g 2 a 2 around each time t n the curvaton crosses around zero, see (7.4). We 
find 

n 4 rr 2 T 2 

g 2 a 2 (t n + St) « - 4 St 2 

! + *(»-§) |t. 
= KUn)5t 2 . (7.7) 



cut 

Thus, the new physical cut-off scale is 

K cut (n)^j - 1 , (7.8) 

closely related to the cut-off scale in the case of m(T) dominated by m a . They are related 
via gyJa^T* = ^gT*/m a A, where A = ^Jgm a a* (see (4.6)). The ratio gT*/m a controls 
how similar the two scales are. 

Following a similar analysis to the case of m{T) ~ m a , we conclude that non-perturbative 
particle production takes place when 

KL(n) > 7rg 2 T 2 (n) (7.9) 
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or equivalently, if 

g > irg T (^) . (7.10) 

This is completely different to the case m(T) « m a , because it does not depend on the 
number of oscillations. Thus, for large g the resonance is initially unblocked and decay 
happens fast, so there is no time to build up the curvature perturbation. However, for 
smaller g, if the resonance was initially blocked, then it remains blocked provided that 
m(T) ~ gT. When the curvaton's thermal corrections become subdominant (T < 7t ran ), 
the calculations in sections 5 and 6 become valid. In that case, we match the solutions 
for a(t) ((3.4) and (7.2)) at the moment when T = T trans , as explained above. This is 
equivalent to making the following substitutions in the equations of sections 5 and 6: 

(TJX \ ( Tfl \ 2 

-yJ(t*, H *^[^r) H *> t^(t-t tran ). (7.ii) 

The non-perturbative effects could then begin either immediately after t = t trans , or they 
could be still blocked for some time. Thus, it is important to determine (a) whether 
gT* > rriu, (b) if this is true, then whether the resonance is blocked, and (c) how decay 
occurs after tt ra n in the regime m{T) ~ m a . 

Note that the above discussion assumed broad resonance, and is therefore only valid if 

«. = ^~^f^f^>l. (7.12) 



So what about narrow resonance in the regime m(T) ~ gT? The narrow resonance becomes 
unblocked only if 



g 2 T 2 (n) + ±qm\Ty cos 2 %^ ^l + 2H*t-l < oj 2 . (7.13) 



Because q < 1 for narrow resonance, this would require g 2 T 2 (n) < g 2 T 2 (n). In other 
words, independent of the number n of curvaton zero-crossings, narrow resonance in this 
regime requires g > g T « 0.3. In section 8, we present the parameter space only for g < 0.1. 

In summary, there are two possible scenarios when the curvaton's thermal mass is 
important. For large couplings, the non-perturbative effects begin immediately and thus 
the correct curvature perturbation is never obtained. For smaller couplings, the non- 
perturbative effects are blocked at least until the temperature drops below Tt rans = m a /g, 
and from there on the analysis for m(T) = m a applies. 



8 Consequences for the curvaton model 

In the above sections we presented analytical results for the curvaton's non-perturbative 
decay in both broad and narrow resonance, for both matter-dominated (MD) and radiation- 
dominated (RD) scenarios. In this section, we first explore how the type of decay varies 
with the parameters g and m a . We then discuss the calculation of the observable C in the 
cases of non-perturbative decay (with and without thermal corrections) and the case of 
perturbative decay. This serves to illustrate the substantial differences in dynamics when 
the realistic mechanism of decay presented in this paper is considered. 
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Figure 2. Showing the type of decay for fixed c* = 10 16 GeV and = 10 12 GeV. For large m CT 
and g, narrow resonance completes in either RD or MD eras (orange and red). For smaller g or 
smaller m a , narrow resonance begins, but does not complete before EWSB (lilac and blue). Above 
the dotted line, the curvaton's thermal corrections can substantially affect the equations of motion. 
Above g = 0.1, the effective mass of the higgs due to interactions with the curvaton is larger than 
the background temperature go* > T*, and thus our formalism for the thermal corrections is not 
applicable. 



8.1 Broad versus narrow resonance 

In order to illustrate the dependence on the parameters of the potential, g and m a , we fix 
cr* = 10 16 GeV and = 10 12 GeV. In this example, narrow resonance always occurs, 
and Fig. (2) shows whether this takes place in MD or RD. Above the dotted line, the 
thermal corrections to the curvaton are initially dominant i.e. gT* > m a . For g > 0.1, the 
effective mass of the higgs is too large for our formalism to be valid. This is because the 
effective higgs mass due to the curvaton-higgs interactions is large, go* > T*. The narrow 
resonance is only efficient for large m a and large g (orange and red). In this case, the 
resonance ends either in MD (smaller g, red) or RD (larger g, orange ). For g < 10" 10 and 
for m a < 10 5 GeV, the narrow resonance is not complete by EWSB (lilac and blue). After 
EWSB, the higgs non-linearities and the term —Xv 2 in the higgs equation of motion may 
no longer be ignored. Three-point interactions are also expected to become important. We 
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leave this calculation for future work. 

In this example, there is no region where broad resonance is important. This is because 
of the thermal corrections to the curvaton, which are largest for large g. These corrections 
cause the decay to be delayed until at least T = m a /g. By the time that these corrections 
become negligible (i.e. T <C m a /g), the resonance parameter q is generally small, and we 
are in the narrow resonance regime. 

As if* decreases (with fixed cr* ) , the region in which our calculations are valid decreases 
slightly. Other than that, the parameter space does not change qualitatively. Changing cr* 
(fixed if*) has a bigger effect. When cr* is small (cr* « if*) then most of the parameter 
space has narrow resonance delayed until EWSB (RD). There is complete narrow resonance 
only if both g and m a are large: g > 1CP 7 and m a > 10 8 . As cr* is increased to Mp, the 
size of the region where narrow resonance completes increases substantially. There is also 
substantially more decay in the MD era. 

Thus, we conclude that narrow resonance dominates the parameter space, and that 
the timescale for this strongly depends on m a , g and cr* , but only weakly on if* . We leave 
a full analysis of the parameter space for future work. 

8.2 The curvature perturbation 

As is well known, the curvature perturbation in the curvaton model is given approximately 
by [l] 

^ _ ff* r dec -g 
37TCT* 



where 

3/9, 

Tdec ~ 



3/V + Ap rad 



(8.2) 

3 decay 



where jdecay = i B R + A? f° r broad resonance and jdecay = j NR + for narrow resonance. 
We consider which regions of the (m a ,g) parameter space can give a sufficiently large value 
of ( to match the observation of £ ~ 10 -5 . We do this by varying cr* and searching for 
the maximum value of £. Typically, £ is small for both small cr* and for large cr*, reaching 
a maximum at some intermediate value. Thus, if the maximum possible £ is larger than 
observations, then it is possible to choose one or two values of cr* that give £ matching 
exactly with observations. However, if this maximum £ is less than the observed value, 
then that point in the (m a ,g, -H*) parameter space is ruled out because there is no choice 
of cr* that gives ( = 10 -5 . 

If we were to consider the curvaton decay width as a free parameter, the maximum £ 
occurs at the border between RD (r^ ec 1) and MD (r^ ec ~ 1) regimes (see (8.1)). In this 
case, roughly estimating this maximum £ gives 

if* 

.V/,,-"" 



Cmax ~ 0.1 Jdecay' 



Observations constrain if* < 10 13 GeV, thus ( max < lCT 6 j 1/4 - This means that jdecay ^> 
10 would be necessary to match observations. In our case, jdecay depends on the parame- 
ters of the model in a highly non-trivial manner. Thus, the estimate (8.3) cannot be relied 
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upon. It does however indicate that very small values of jdecay are not able to produce 
sufficient £. This can be understood because the curvaton must be subdominant during 
inflation, and then requires a substantial period of growth to produce a sufficiently large 
curvature perturbation. 

In the scenario considered in this paper, we emphasise that non-perturbative decay 
in the presence of thermal corrections is the only physically relevant mechanism for de- 
cay. However, to illustrate the substantial difference in the parameter space when this 
non-perturbative decay is correctly accounted for, we present calculations of £ in three 
scenarios: 



(i) Perturbative decay. This can be parameterised by an effective decay width as 

where A is an effective coupling. In this case, the maximum £ is 

H 1 

Ua^O-lj^-^- (8.5) 

Thus, a perturbatively decaying curvaton requires A < 10 _4 (-£/*/10 12 GeV) 2 to match the 
observed £. Note that it is never correct to ignore the possibility of non-perturbative decay, 
and we only give this result for comparison with earlier work. 

(ii) Non-perturbative decay in the absence of thermal couplings. As shown in 
section 4, this has jdecay — C(l — 100) and in this case it is never possible to obtain a large- 
enough (. Thus, an important conclusion is that in the absence of a thermal background, 
a curvaton model with this type of non-perturbative decay is ruled out. This scenario ap- 
pears contrived, however it would in principle be applicable if the inflaton decayed to some 
hidden sector and then the curvaton became dominant, before decaying into SM particles. 

(iii) Non-perturbative decay including thermal couplings. In this case the 
effective decay width is calculated using results from sections 5, 6 and 7. As we have 
discussed, we do not consider the process of decay after EWSB. In that case, it would be 
important to check whether decay occurs early enough to avoid spoiling the predictions of 
big bang nucleosynthesis. If we assume that decay happens instantly at EWSB (unlikely 
to be a good assumption), then we have checked that we can obtain £ = 10~ 5 in almost 
the entire parameter space presented in Fig. (2), by tuning ex* to suitable values. For 
if* = 10 12 GeV, only a combination of large m a > 10 10 GeV and very large g > 10~ 5 
is ruled out. Thus, in contrast to the case (ii) above, a viable model is possible even for 
relatively large values of the curvaton-higgs coupling. This important conclusion applies 
to any curvaton model where the decay products are present in the thermal background. 

As decreases, the parameter space also decreases, until for ~ 10 9 GeV, there 
appears to be no parameter space remaining for g > 10~ 15 . This decrease is common to 
all curvaton models and it occurs because £ oc (see (8.1)). The exact parameter space 
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depends on the details of the decay after EWSB, which we do not discuss in this paper. 
Also note that non-Gaussianity is likely to be very different in cases (ii) and (hi), even if 
both scenarios can produce the correct £• This should be investigated numerically 

9 Discussion and conclusions 

We have investigated a simple, realistic and concrete model for the decay of the curvaton 
condensate. In the model the curvaton a is coupled to the SM Higgs field through a 
(j 2 <j)t$ term, which is the only renormalizable coupling of a singlet — like the curvaton 
- to the SM. With such a coupling, there is no direct, tree-level decay channel available. 
Instead, curvaton decay takes place by resonant production of higgs particles. This is a 
non-perturbative effect, which typically is highly efficient and could make the curvaton 
decay very quickly, thereby preventing the initially subdominant field perturbation from 
growing and giving rise to the observed curvature perturbation when the curvaton conden- 
sate decays. However, here we show that the thermal background due to the inflaton decay 
products completely changes the dynamics of the non-perturbative production of higgs 
particles. This is in contrast to studies of inflationary (p)reheating [7-11], where there is 
no pre-existing thermal background. The effect of a thermal background on curvaton decay 
has recently been addressed in a different context in [16] . 

We assume that when the curvaton oscillations start, the inflaton has decayed into SM 
degrees of freedom, which have thermalised among themselves. This generates an effective 
thermal mass both for the curvaton and for the higgs. As a consequence, as we show in this 
paper, resonant production of higgs particles is blocked by the thermal higgs mass, which 
is proportional to the temperature of the thermal bath. Under these circumstances, the 
eventual onset of the resonant higgs production depends on various effects: on the nature 
of the resonance; on the expansion law of the universe; and of course, on the curvaton 
mass m a , the curvaton- higgs coupling strength g, and the value of the inflationary Hubble 
rate H*. The effective decay rate of the curvaton condensate, or the actual time when 
a given fraction of the initial condensate energy has been converted to SM particles, also 
depends on the effectiveness of the energy transfer once the non-perturbative effects become 
unblocked. 

All of these issues were discussed in detail in sections 5, 6 and 7. There we treated 
separately the cases of broad and narrow resonances, noticing that the curvaton usually 
starts its oscillations in the regime of broad resonance, which then can evolve into a narrow 
resonance while resonant higgs production remains blocked. We also discussed separately 
the radiation dominated universe and the matter (curvaton-oscillations) dominated uni- 
verse. We found that for a range of parameters, the curvaton decay can be delayed all the 
way down to the electroweak symmetry breaking scale. Our main results are the equations 
(5.12), (5.18), (6.8) and (6.7) which define the time scales at which non-perturbative effects 
first become important, and equations (5.24) and (6.12), which describe the time scales 
for efficient decay after resonant production has begun. We also discussed the thermal 
correction to the curvaton's mass. This affects the curvaton and higgs dynamics in a sub- 
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stantial fraction of the parameter space, and tends to add a further delay to the decay of 
the curvaton condensate. The main equations in this case are (7.8), (7.10) and (7.11). 

Compared to the usual inflaton (p)reheating, the existence of a thermal background sig- 
nificantly impacts the whole non-perturbative particle production process. A similar situa- 
tion could be encountered in other cosmological scenarios involving the SM degrees of free- 
dom as decay products, such as (p)reheating after Higgs Inflation [17-19] or (p)reheating 
after MSSM flat-direction inflation [20-22]. In both of these models, the inflaton decays 
by non-perturbative effects into SM or MSSM degrees of freedom. Because these degrees 
of freedom are expected to thermalise very quickly, the results presented in this paper 
could have an important impact on their dynamics. In addition, our results could have 
a direct impact on gravitational wave backgrounds created during (p)reheating [23-27]. 
The delay of the resonance due to thermal effects means that the peak frequency of the 
gravitational wave backgrounds would be shifted towards much smaller values, favouring 
their observability. 

In section 8 we considered the type of decay occurring in the space of the curvaton 
parameters, for fixed H* and fixed <7* — the result is displayed in Fig. (2). Decay by 
narrow resonance dominates the parameter space, and occurs either in matter domination 
or in radiation domination. For small g or small m a , narrow resonance begins but does 
not complete before EWSB. Under various assumptions, including assumptions about the 
decay after EWSB, we checked that the observed £ = 10 -5 can be obtained in most of 
the parameter space, except for a region with very large m a and large g. We find that 
realistic couplings g < 10 _1 give a viable curvaton model. This is in contrast to the case 
of non-perturbative decay without a thermal background. In that case, £ = 10 -5 is not 
possible because the curvaton decays too quickly, and the model is ruled out. 

In the case where relatively large values of g are allowed, the Colemann- Weinberg 
radiative corrections to the curvaton potential may become important and have a significant 
impact on the results [28, 29]. Here we only considered the simplest quadratic curvaton 
potential. Future work could consider both the addition of a quartic term Act 4 and a 
Coleman- Weinberg term. In future work, it will also be necessary to calculate the non- 
Gaussianity produced in this scenario. 

In summary, we have carefully calculated the time scales involved in the decay of 
the curvaton via resonant production of Higgs bosons. In particular, we showed that the 
thermal corrections to the higgs mass substantially delay the decay of the curvaton. This 
means that 'natural' values of the curvaton-higgs coupling can give rise to the observed 
perturbation amplitude without fine-tuning of the decay rate. 

Note added: After completing this work, we were made aware of [30], which addresses 
similar issues in a different context. 
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A Timescale of energy transfer in broad resonance 

In broad resonance, the timescale of energy transfer is short enough to be effectively negligi- 
ble in our calculations, that is Aj <C j BR . In this appendix, we demonstrate the calculation 
of the timescale for both in RD and in MD, in the case when the curvaton effective mass 
is m(T) ~ m a . The calculation is substantially different compared to the usual inflaton 
(p)reheating case, because of the background thermal ensemble. 

We now switch to a more compact notation for the j-dependence of quantities. In the 
new notation, nfcC/ B R + Aj) becomes njj. , a (j BR + Aj) becomes a(Aj)> an d similarly for 
other quantities. Let us recall (5.22), which in the new notation is given by 



n (Ai) 
k 



H(1 + 2W^) 



< + (i + 2< vi o) ) + e ( n a + 2w ^) < 



b=2 \c=b+l 



k ' 

(A.l) 



which describes the higgs occupation number Aj semi-oscillations of the curvaton since the 
broad resonance became unblocked. The second term (containing index b and c) on the 
right hand side of (A.l) describes the accumulated particles created out of vacuum at each 
curvaton zero-crossing, as well as the stimulated creation of particles out of the previously 
created particles from vacuum. This term becomes subdominant after few oscillations, and 
thus will be neglected. Within the remaining term (containing index a), the contributions 
to the occupation number come from: (i) the first creation of particles out of vacuum, de- 
scribed by wf^\ and (ii) the deformation of the original thermal higgs ensemble, described 
by (1 + 2W^)nf \ where 

n i°^ = nt k(T{j w )) is the thermal ensemble occupation number 
at the moment when resonant production begins. 

A.l Deriving the characteristic time-scale Aj for efficient energy transfer 

To find the energy transferred into the higgs from the curvaton, we should solve (A.l) 
recursively to obtain the total energy density of the higgs ensemble pjj(Aj), given by 

Ph (Aj) ee J dktfnMy/kyafa+ml^, (A.2) 

where Tn>h,(Aj) if the higgs effective mass at j = j BR + Aj (which we discuss later). However, 
there is a simple approximation that by way of comparison of characteristic scales effectively 
removes the need for this recursive solution, which we now describe. First, note that the 
spectral shape of Wu~ is different from that of nl . The former is the window function 
given by (5.7), 

W W = e -<^)(%)/e) 2 e -( fc / fe -) 2 , (A.3) 
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and the latter is a Bose-Einstein distribution, 



11 



(o) 



l k2 / a (i)+ m t(i)/ T W 



(A.4) 



In the infrared (IR) limit (k <C k^ t ), both distributions are approximately constant. In 
the case of the thermal distribution this is because, as we demonstrate in section A. 2, the 
higgs mass is much larger than the physical cut off, rn 2 h ^ S> (K^ t ) 2 . In the ultraviolet 
(UV) limit, both distributions are exponentially suppressed 5 . However, this occurs at well 



separated scales: W\ is suppressed for k > k c J t , whereas (1 + 2W k 1 )n k is suppressed for 
^ a (i) m h,(i) ^ k^cut (t ne calculation of m 2 h is shown in Sec. A. 2). However, although 



a (i) m h.(i) ^ ^cut' the term 11^2 (■"■ + ) wri ich multiplies both distributions is also 
suppressed for k > k^ t . Thus, there is only one characteristic scale, given by k^ t . The 
occupation number is then 



rAj 



n 



(Aj) 



Aj 



ri( i+2 ^ 



WP + (1 + 2W, 



(!)^(°) 



n, 



oc const., 



for k < fcW , 



(A.5) 



whereas for k > k^ t is exponentially suppressed. The occupation number can then be 
approximated as 



n 



(Aj) 



/(<?) 



Aj 



U( 1+2W k 



w. 



(1) 



where 



/(<?) = 1 + 



(1 + 2W, 



n 



W, 



(i) 



(A.6) 



(A.7) 



is defined only for k <C fc^. Thus the function f(q) quantifies the difference in amplitude 
between the vacuum- and thermal-originated occupation number distributions. It is scale 
independent and only depends on the resonant parameter q via the mass term in the 
thermal distribution ni . We can obtain f(q) by evaluating the expression above at any k 



smaller than k c J t . Plugging in the expressions of W k and n k and, choosing k sufficiently 
small as compared to k^ t , we can explicitly write 

f(q) ra 1 + 



e 9 T q 



e + 1) 

1/4 



(Ai 



» (k 



cat 



which we now demonstrate. The 



A. 2 Calculation of m h ^j) 

Our reasoning above assumed cJ(Aj) m h, (Aj) 
effective higgs mass contains contributions both from the interactions with the curvaton 
and the interactions with the thermal bath. Averaging over a curvaton oscillation period 
we obtain 

m 2 h (Aj) = -g 2 Z 2 +g 2 T 2 (A.9) 



3 With different powers of k. 
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where ^(Aj) is given by (3.7) in RD or (3.15) in MD. The blocking is lifted (i.e. resonant 
production begins) when g^a^^^T^^/k^p) 2 < 7r _1 ~ O(0. 1). The ratio between the 

-y* 

mftAj) 



higgs mass and the typical physical momentum = k^t / a (Aj) i s 



(€P/a {Aj) ) 2 



cut 



2 n A V z 

g z a (Aj)^(Aj) 

2 {k (^h2 



2m CT 

aV 2 



+ 0(0.1) 



+ 0(0.1) 



(A.10) 



Note that this relation is true in both RD and MD backgrounds. 



We have thus demonstrated that CL 2 A ^m 2 i 
in broad resonance, a ^v/) m h(Aj) (^cuP) 2 

,2 U/2 



( ? 1 ^ 2 (^cu/^) 2 i n general and, then 



factor (k /a 2 A jj + m 



S> {k^i 1 ) 2 , as we claimed. We see from this that the 
in (A. 2) is totally dominated by the mass contribution, as 



it correspond to non-relativistic higgs particles. That is, broad resonance always gives rise 

v cut 



to non-relativistic decay products. This allow us to write Jk 2 + af 



(Aj) m h,(Aj) 



in (A. 2). 

A. 3 Energy transfer 

We are now ready to compute the energy transfer into the higgs in the case of broad 
resonance and in the presence of a thermal ensemble. This is quite similar to the calculation 
in section 4. 



Due to the condition of non-thermal blocking g 2 a 2 ^T 2 ^ 



W, 



w = e -^4 o w( T w/*St) a e -^(*/fcS*) 



kcnt/v, then 

-<*/e) 2 



(A.ll) 



We Taylor expand the window functions Wl, following a similar approach as in the 
case of no thermal corrections (section 4.2), and obtain 



n 



(Aj) 



H(i+2w, 
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(A.12) 
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with 

where we have used k^ t /k^ t « 1 because it only depends very weakly on j. 

Notice that (A. 12) has all the properties that we previously discussed: it is constant 
for k <C k^ t and it is exponentially suppressed for k S> k^ t . 

We can finally substitute (A. 12) into (A. 2), and use the higgs mass definition from 
section A. 2 to obtain 

„<4fl„4/(«) {, , 2^- 1 (A£}) 4 
Ph 



e 



1+ 2 ) / dxx 2 e-^ A ^y/x 2 + q 1 / 2 (A.14) 



with the factor 4 accounting for all the higgs components (which follow the same particle 
creation dynamics) and where we have defined x = k/k^ t . Because the integrand is 
exponentially suppressed for x > 1, we can approximate sj x 2 + q 1 / 2 m q l l A , as already 
mentioned before. The integral can then be easily estimated analytically, and we find 

4 /(g) ^ , 2^- 1 (*2)* g 1 / 4 

4, 



eF^(Aj) \ 1+ e) (2tt) 3 ' (A-15) 



where we have used 47r J dxx 2 e nx2 = 1. This is the expression in (5.25). 
B Timescale of energy transfer in narrow resonance 

The mechanism of energy transfer in a narrow resonance is different from the case of broad 
resonance, as we now discuss. We consider the case where the curvaton effective mass is 
m(T) ~ m a . 

As commented in the main text, we can ignore the higgs non-linearities provided the 
energy transfer occurs before EWSB. This is our basic assumption for this calculation. The 
higgs component mode equation, written as in section 6 with Xa = a 3 ^ 2 4>a, is given by 

Xa + {^2 + 9 2 T T 2 + 4q sin 2 (m a t + tt/8)^ X a = . (B.l) 

Here q = g 2 T, 2 /Aml < 1. 
B.l Mathieu equation 

Defining a dimensionless time variable z = (m^t + vr/8), we can re- write (B.l) as a time- 
dependent "Mathieu equation" : 

d 2 Xc 



dz 2 



+ [A k (t)-2q(t)cos(2z)] Xa = 0, (B.2) 
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with Ak(t) = ^0j^v — I- 9^(T(t)/m a ) 2 + 2q(t). The solutions to the time-independent 
Mathieu equation are well known. Because there is no violation of adiabaticity for q < 1, 
we can consider an 'instantaneous' Mathieu equation at each time, and thus we can use 
the known band structure of the Mathieu equation in the narrow resonance regime. 

Solutions to the Mathieu equation are oscillatory and either stable, or are have ex- 

'.(').' 



ponentially unstable amplitudes growing as Xa{k,t) oc exp y[r k J z\. The growing solutions 
occur in a set of narrow resonance bands labelled by I. Within these momentum bands, 
nk{t) oc exp (2[a^z\. The narrow resonance bands are characterised by ~ I 2 , with 



I = 1,2,3, .. .. Each band in momentum space has a width dictated by AA k ~ q l . The 
first band I = 1 is the widest and most important one. For this band, Ak ~ 1 ± q and 



{k ^+9 2 T (T/ ma f = l-2 q ±q. (B.3) 



0/ 



Non-perturbative resonant production will begin below the temperature T/vb> when 
(B.3) has non- negative solutions for k 2 . Thus, choosing the central value of the band, 
i.e. taking (1 — 2q) on the right hand side of (B.3), and setting k = on the left hand side, 
gives 

T NR ^ { -^m a , (B.4) 
which coincides with (6.6). The width of the first band is 

— ~ q 1 / 2 m a . (B.5) 
a 

Thus, only very IR modes with k < Ak will be excited. It is interesting to note that 
without thermal corrections, a set of modes with larger momentum would be excited, with 
k/a ~ Trio- (1 -?±(?/2). The interactions with the thermal bath reduce the allowed phase 
space and make the production of particles inefficient, only allowing for very IR modes to 
be produced. This reduced phase space means that the relevant initial thermal distribution 
has a simple form. The initial spectrum of the occupation number nP , should be the Bose- 
Einstein distribution. However, applying the band condition (B.3) and recalling that the 
effective higgs mass is rr? h = g 2 T 2 + ^-S 2 ~ 2qm 2 , we find 

-i 



(0) = ^(k/a^+g^+2qml/T _ l 

= ^ e 9 T ^(fc/m CT ) 2 /a 2 +4(T/m CT )2+2q _ ^ ^ 



+ (B.6) 

In other words, the initial condition for the modes within the resonance band is a temperature- 
independent distribution with scale-invariant amplitude rcj? ~ (c 9t — l)" 1 ~ 3. 



- 35 - 



B.2 Energy transfer 

The first resonance band has Floquet index given by [10] 



( \(k/a) 2 + g 2 T 2 ] 1/2 

t /2) { a-2,Xl (R7) 



The maximum of this is in the centre of the band ((k / a) 2 + g 2 T 2 = (1 — 2q)m 2 ), and gives 
/x™ ax = q/2. The Floquet index very quickly goes to zero as k approaches the band width, 
Ak = aq 1 l 2 m u . Thus, we approximate the Floquet index as 

f % , k < Ak/2 
Mfc ~ { ~ (B.8) 

\ , k > Ak/2 V ' 

This gross approximation is sufficient to capture the essence of the physics: only very IR 
modes will be excited and the typical amplitude of the Floquet index is of order ~ q/2. 

We are then ready to estimate the energy density of the Higgs field once the narrow 
resonance process begins to take place. Counting the number of zero crossings of the 
curvaton asz = irAj since the end of the thermal blocking, the excited modes grow out of 
the initial Bose-Einstein distribution as 

Including again a factor 4 to account for all higgs components, we then obtain 

= I d kk 2 w - Y) ^k,a) 2 + g 2 T T 2 + 2qrnl 

2m a e *«4j f Ak / 2 



dkk 2 ^l±q 



i pi"gAj 

qV 2 {l + 0(q))mi (B.10) 



3(2vr) 2 (e fl T - 1 

which is the formula in (6.11) that we provided in the core of the text. 
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